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ABSTRACT 


The  purpose  of  this  first  phase  of  a  study  on  the  morphology  of  planetary 
waves  was  to  estimate  global  temperatures  for  the  upper  stratosphere  and  lower 
mesosphere  from  the  radiances  of  the  SCR-B  experiment  on  Nimbus  5.  This  has 
been  successfully  achieved  using  a  regression  of  the  radiances  on  meteorological 
rocket  data.  In  general,  the  accuracy  of  the  retrieved  temperatures,  i.e. 
standard  deviation  of  the  error,  varies  from  about  3°K  at  10  mb  to  about 
5.5°K  at  0.4  mb  as  compared  to  meteorological  rocket  values.  Layer  mean  temper¬ 
atures  are  about  one-half  degree  more  accurate  than  temperatures  at  the  border¬ 
ing  levels.  Details  of  the  procedure,  which  included  extensive  prior  prepara¬ 
tion  of  the  radiances,  are  given  here.  Global  temperatures  at  10,  5,  2,  1 
and  0.4  mb  along  the  orbit,  and  the  mean  temperatures  (thicknesses)  between 
the  levels,  are  now  ready  to  apply  to  the  next  phase  of  this  study,  which  will 
obtain  heights  of  those  surfaces  at  lat i tude- longi tude  gridpoints.  The  planetary 
waves  will  be  computed  from  these  gridpoint  values. 
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INTRODUCTION 


TI»?  past  several  years  have  seen  a  tremendous  surge  of  interest  in  the 
upper  atmospheric  regions  above  10  mb.  Better  knowledge  of  the  dynamics  and 
thermodynamics  of  these  regions  is  now  seen  to  be  essential  to  improved  under¬ 
standing  of,  for  example,  ozone  photochemistry  and  flux,  the  morphology  of 
stratospheric  planetary  waves  and  warmings,  and  possible  links  between  upper 
and  lower  atmospheric  phenomena.  In  particular,  temperature  perturbations  in 
the  region  may  be  useful  precursors  of  perturbations  in  the  troposphere,  i.e. 
weather  anomalies. 

One  of  the  first  requi rements  of  such  research  is  a  daily  data  set  which 
is  much  more  complete  in  spatial  coverage  than  any  previously  available  data 
set.  The  purpose  of  this  first  phase  of  the  study  was  to  provide  global 
temperatures  above  10  mb  from  the  SCR-B  radiometer  on  Nimbus  5  for  1973-1974. 
In  view  of  the  experience  and  improvements  gained  in  processing  SCR-B,  the 
SCR-A  data,  from  Nimbus  4,  1970-1972,  have  been  processed  in  the  same  manner 
so  that  the  entire  data  set  is  consistent.  A  description  of  the  processing 
used  in  both  sets  is  included  here  for  convenient  reference.  Together,  the 
SCR-A  and  SCR-B  data  now  provide  a  nearly  5-year  data  set  from  which  many 
studies  can  be  made. 

IX.  SELECTIVE  CHOrPER  RADIOMETER  DATA 

A .  The  SCR  Instrument 

From  April  1970  through  December  1974,  two  instruments,  SCR-A  and 
SCR-B,  were  flown  on  NASA's  polar-orbiting  satellites.  Nimbus  4  and  Nimbus  5, 
respectively.  The  SCR  (Selective  Chopper  Radiometer)  instruments,  which  are 
described  fully  in  References  1  and  2,  were  designed,  in  part,  to  sense  radia¬ 
tion  upwelling  from  regions  10  mb  to  0.4  mb.  The  weighting  functions  for  both 
instruments  are  shown  in  Figure  1;  the  radiances  which  are  available  for  de¬ 
termining  temperature  in  the  regions  above  10  mb  are  from  Channels  A  and  B 
in  SCR-A,  and  Channels  B12,  B23,  B34,  and  B4  in  SCR-B.  These  radiances  have 
been  "cleaned"  carefully,  regressed  with  coincident  rocketsonde  data,  and 
used  to  produce  a  56-month  time  series  of  sub-orbital  temperatures  at  10,  5, 

2,  1,  and  .4  mb  and  mean  temperatures  (thicknesses)  in  the  layers  between 
those  levels. 


B.  Data  Sources 


Most  of  the  SCR-A  and  SCR-B  data  tapes  were  received  from  the  National 
Space  Science  Data  Center  (NSSDC).  These  were  copies  of  tapes  supplied  to 
NSSDC  by  the  Oxford  experimenters.  The  48  SCR-A  tapes  covered  the  period  27 
July  1970  to  30  January  1973.  The  period  11  April  1970  to  27  July  1970  was 
covered  by  a  set  of  radiance  tapes  available  only  at  CDC.  These  data  were 
calibrated  at  CDC  some  time  ago  from  the  SCR-A  raw  radiances  extracted  from  a 
set  of  800  Sensory  Data  Tapes  supplied  by  NASA.  The  63  SCR  B  tapes  covered 
the  period  13  December  1972  to  26  December  1974.  Formats  of  the  NSSDC-supplied 
tapes  are  given  in  References  3  and  4.  The  two  sets  of  data  tapes  were  organized 
similarly,  but  the  details  of  the  formats  were  quite  different  as  the  instru¬ 
ments  themselves  were  quite  different. 

C .  SCR-B  Radiances 

Examples  of  the  calibrated  radiances  from  the  SCR-B  data  tapes  are 
shown  in  Figures  2  and  3.  The  two  orbits  plotted  are  nearly  complete  from  the 
Northern  Hemisphere  (NH)  descending  into  the  Southern  Hemisphere  (SH)  and  then 
ascending  into  the  NH  again.  The  data  were  consistent  enough  from  point  to 
point  so  that  it  was  necessary  to  plot  only  every  fifth  point.  The  radiances 
from  the  B-difference  channels  (B13,  B23,  B34)  were  previously  smoothed  by 
the  Oxford  experimenters.  They  used  a  five-point  running  mean  and  applied 
it  to  the  orbital  data  four  times.  Details  of  the  calibration  processing 
are  given  in  Reference  2. 

D.  SCR-A  Radiances 

Figures  4  and  5  show  radiances  from  the  two  highest  channels  of 
SCR-A  at  every  point  along  one  orbit  in  December  1971.  The  most  striking 
feature  of  the  SCR-A  radiances  is  the  high  frequency  noise,  particularly  on 
Channel  A.  These  radiances  were  not  smoothed  during  the  calibration  proces¬ 
sing  at  Oxford  as  the  SCR-B  radiances  were.  Channel  A  and  Channel  B  were 
therefore  smoothed  here  with  the  filter  which  is  discussed  later. 

III.  DATA  EXTRACTION  AND  SORT 

Preliminary  processing  of  the  SCR  data  from  the  two  satellites  was  done 
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separately,  with  two  sets  of  programs.  The  input  data  tape  formats  were  differ¬ 
ent  primarily  because  SCR-B  had  more  and  different  channels  than  SCR-A.  In 
addition,  the  SCR-A  radiances  had  to  be  filtered  along  the  orbits.  However, 
most  of  the  techniques  described  below  were  used  on  both  SCR-A  and  SCR-B  data. 

A.  Extracting  and  First  Sort 

The  calibrated  radiances  were  just  a  small  part  of  the  data  on  the 
tapes  received  from  NSSDC.  The  tapes  also  contained  identification  words, 
engineering  data  concerning  the  state  of  the  instrument  and  the  satellite, 
calibration  parameters,  and  data  from  other  experiments.  The  data  extracted 
and  saved  on  a  new  set  of  tapes  included  orbit  number,  day,  time,  latitude, 
longitude,  radiances,  and  24  flag  bits.  Some  of  the  problems  encountered  in¬ 
cluded  sometimes  unreliable  orbit  number,  day  number,  or  time.  In  some  cases, 
the  latitude  and/or  longitude  were  obviously  incorrect.  Sometimes  certain 
identification  records  were  either  missing  or  duplicated,  making  it  difficult 
to  distinguish  ends  of  orbits  and  ends  of  days.  The  biggest  problem  was  the 
data  being  out  of  order.  Before  the  data  could  be  filtered  they  had  to  be 
arranged  chronologically.  The  first  re-sort  was  done  by  copying  the  extract 
tapes  and  re-ordering  individual  orbits  or  whole  days  that  were  obviously  out 
of  order. 

B.  Merge  CPC  Calibrated  Radiances  for  SCR-A 

A  list  of  the  data  missing  from  the  re-sorted  SCR-A  tapes  was  compared 
with  an  inventory  of  the  data  that  had  been  sorted  and  calibrated  at  CDC.  The 
data  available  at  CDC  for  filling  in  gaps  in  the  NSSDC  data  included  every  other 
day  for  the  first  three  months  of  SCR-A  and  several  days  scattered  through  the 
first  year.  The  new  data  were  written  on  another  tape  in  the  same  format  as 
the  data  extracted  from  the  NSSDC  tapes,  with  a  flag  to  denote  the  source  (NSSDC 
or  CDC).  In  addition,  the  CDC-cal ibrated  data  had  to  be  adjusted  to  reflect 
some  changes  in  the  calibration  procedure  which  were  discovered  after  the  orig¬ 
inal  calibration  had  been  completed.  The  extracted  data  from  the  two  sources 
were  then  merged  by  copying  to  a  new  set  of  tapes.  This  step  also  took  care 
of  a  few  orbits  that  had  been  missed  in  the  previous  re-sort. 

C.  Final  Sort 


The  final  sorting  technique  was  designed  to  handle  out-of-order  and 
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duplicate  data  within  orbits.  It  also  handled  one  or  two  orbits  near  midnight 
GMT  that  were  almost  always  out  of  order  and  that  were  not  re-sorted  previously. 
For  each  24-hour  period  a  5400  word  array  was  generated.  There  was  one  slot 
in  the  array  for  each  16-second  satellite  observation  interval,  called  a  major 
frame.  If  the  time  corresponding  to  the  first  slot  is  known,  the  times  of  all 
slots  are  known  since  each  one  is  16  seconds  after  the  one  before.  The  array 
was  positioned  to  cover  the  period  from  12Z  one  day  to  12Z  the  following  day. 
Since  the  out-of-order  orbits  near  00Z  were  then  near  the  center  of  the  array, 
it  was  not  necessary  to  switch  back  and  forth  between  arrays  for  different  days. 

The  5400  words  were  first  cleared  to  zero  which  represents  "missing". 
Words  corresponding  to  major  frames  without  good  data  remain  zero.  Each  major 
frame  was  then  examined  to  determine  if  there  were  data  worth  keeping.  A  major 
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frame  was  rejected  if  any  of  the  radiances  were  missing,  or  if  the  time  or  posi¬ 
tion  was  obviously  bad  (e.g.  hour  was  greater  than  23,  or  latitude  was  greater 
than  80°N,  or  longitude  was  greater  than  360°E) .  The  following  parameters  were 
saved  from  a  good  SCR-A  major  frame:  latitude,  longitude,  Channel  A  radiance, 
Channel  B  radiance,  and  a  flag  to  distinguish  the  source.  The  SCR-B  parameters 
saved  were  latitude,  longitude  and  radiances  for  Channels  B12,  B23,  B34,  and 
B4.  The  position  of  the  major  frame  in  the  5400  word  array  was  determined  from 
the  time,  and  the  parameters  were  stored  in  the  appropriate  word.  When  a  new 
day  was  encountered,  the  5400  word  array  was  written  to  tape  and  the  process 
was  repeated  for  the  new  day.  The  generated  tapes  were  very  compact,  having 
two  year's  data  per  reel. 

IV.  ERROR  CHECKING  AND  FILTERING 
A .  SCR-B  Error  Checking 

The  SCR-B  data  were  examined  to  identify  and  remove  erroneous  values. 
Latitudes  and  longitudes  for  approximately  ten  orbits  in  the  two-year  period 
were  found  to  be  incorrect.  These  incorrect  locations  were  identified  by  compar¬ 
ing  the  ascending  and  descending  nodes  (equator  crossings)  in  the  data  with 
nodes  computed  independently  based  on  the  regularity  of  the  Nimbus-5  orbits. 

No  attempt  was  made  to  correct  the  erroneous  locations  as  the  amount  of  data 
va s  sma 1 1 . 
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The  data  series,  which  were  normal ly  quite  smooth,  had  occasional 
random  spikes,  which  were  identified  and  renoved  by  another  screening  tech¬ 
nique.  Minima  and  maxima  of  all  radiances  in  10°  latitude  belts  were  computed 
for  each  day.  It  was  found  that  extremes  for  each  latitude  belt  are  very  con¬ 
sistent  from  day  to  day  except  when  spikes  in  the  data  caused  an  unusual  minimum 
or  maximum.  By  examining  the  time  series  of  minima  and  maxima  it  was  possible 
to  specify  a  range  of  allowable  radiances  for  each  channel  for  each  day  and 
latitude  belt.  Any  radiances  outside  the  allowable  ranges  were  removed.  This 
technique  typically  removed  fewer  than  five  points  per  day. 

B .  SCR-A  Error  Checking 

The  same  error  checking  procedures  just  described  for  SCR-B  were 
also  applied  to  SCR-A.  A  greater  amount  of  data  was  determined  to  be  bad  than 
was  found  for  SCR-B.  Seventeen  orbits  were  removed  because  of  bad  locations. 

The  technique  of  specifying  ranges  of  allowable  radiances  also  worked  for 
SCR-A,  even  though  the  data  were  much  noisier. 

C .  SCR-A  Filtering 

Following  this  error  screening,  it  was  still  necessary  to  filter 
the  SCR-A  data  to  remove  the  noise  described  in  Section  II-C.  Since  Channels 
A  and  B  were  to  be  used  together,  the  same  filter  was  used  for  both  channels, 
although  A  (Figure  4)  was  obviously  more  roisy  than  B  (Figure  5).  Figure  6 
shows  Channel  A  radiances  from  part  of  one  orbit  before  and  after  filtering. 

A  simple  running  mean  removed  much  of  the  random  noise  evident  in  the  ur.filtered 
data.  It  was  decided  to  use  a  five-point  running  mean,  applied  four  times,  to 
be  consistent  with  the  smoothing  done  at  Oxford  on  the  SCR-B  data.  As  mentioned 
in  the  Oxford  documentation,  this  is  essentially  a  Gaussian  filter  with  a 
half-width  of  about  seven  points  (major  frames).  The  chosen  filter  is  shown 
as  the  heaviest  line  in  Figure  7.  Figure  8  shows  how  effective  this  17-point 
filter  is  in  removing  the  noise  in  Channel  A  radiances.  The  three  short  seg¬ 
ments  shown  are  from  the  orbit  plotted  in  Figure  4. 

Before  the  SCR-B  data  were  smoothed  by  Oxford,  short  gaps  were  filled 
by  linear  interpolation.  Tests  were  made  to  determine  the  effects  this  inter¬ 
polation  would  have  on  the  SCR-A  data.  Five  points  in  the  center  of  each  data 
segment  of  Figure  8  were  removed  and  replaced  by  linear  interpolation.  The 


results  after  filtering  the  interpolated  segments  are  shown  by  the  dots  in 
Figure  9.  The  crosses  show  the  results  after  filling  the  gaps  with  the  filter 
instead  of  by  linear  interpolation,  i.e.  a  missing  point  was  filled  by  the 
weighted  average  of  the  nearby  non-missing  points.  The  solid  curves  are  the 
filtered  results  from  Figure  8  plotted  with  a  twice  expanded  ordinate.  The 
fiLtered  values  at  the  interpolated  points  are  quite  different  from  the  fil¬ 
tered  values  when  the  points  are  not  missing.  Also,  four  points  on  either  side 
of  the  gaps  are  strongly  affected.  The  filtered  values  for  the  gap  and  four 
points  on  each  side  are  quite  different  from  the  original  smoothed  data. 

For  the  above  reasons,  when  the  SCR-A  data  series  was  actually  fil¬ 
tered,  gaps  were  not  interpolated  because  of  the  possibility  of  introducing 
large  errors.  Data  near  gaps  must  be  used  cautiously  because  of  possible 
errors  propagating  from  the  gaps.  Rather  than  throwing  out  the  data  neigh¬ 
boring  the  gaps  during  the  filtering  process,  a  flag  was  added  to  the  data  for 
each  major  frame.  This  flag  showed  the  number  of  points  (up  to  17)  and  the  sum 
of  the  filter  weights  corresponding  to  the  non-missing  points  under  the  filter. 
The  filtered  output  was  also  written  as  a  5400  word  array  for  each  24  hour 
period.  A  few  data  points  were  removed  because  both  Channel  A  and  Channel  B 
were  not  present. 

SCR-A  data  after  17  April  1972  were  filtered  separately  because,  on 
that  date.  Channels  C  and  D  were  lost.  Channel  C  had  been  used  in  the  calibra¬ 
tion  of  Channels  A  and  B  before  that  data,  and,  with  Channel  C  gone,  the  char¬ 
acter  of  the  calibrated  radiances  for  Channels  A  and  B  changed  abruptly. 

V.  REGRESSIONS 

After  the  radiances  had  been  sorted  chronologically  and  filtered,  meteor¬ 
ological  parameters  such  as  temperatures  and  thicknesses  could  be  derived  from 
them.  The  most  reliable  way  to  determine  temperatures  and  thicknesses  from 
radiances  is  by  regression,  as  described  in  Reference  5.  Regression  coefficient 
were  based  upon  a  coincident  set  of  meteorological  rocket  network  (MRN)  and  SCR 
data.  The  13  MRN  stations  used,  primarily  from  the  NH,  are  described  in  Table  1 

Rocket  soundings  were  taken  about  once  a  week  at  each  station,  but  not 
usually  when  the  satellite  was  passing  over  the  station;  therefore,  the  number 


of  coincident  observation:,  was  quite  small.  Tile  sample  size  was  increased  by 
forming  a  time  series  of  all  available  satellite  radiances  at  each  of  the  rocket 
stations,  as  discussed  m  tin-  next  section.  From  the  lime  series  at  a  £ i von 
station,  radiance  values  were  interpolated  in  time  to  coincide  exactly  with 
the  times  of  the  rocket  firings.  These  space/time  "coincident"  data  provided 
the  data  set  for  the  regressions. 

A .  Time  Set  io:.  tieneiat  ion 

Nimbus  4  and  Nimbus  !>  were  local  noon/locai  midnight  satellites.  Any 
point  on  Earth ,  81)' N  to  80' S ,  had  one  or  more  nearby  daytime  overpasses  and 
one  or  more  nearby  nighttime  overpasses.  As  an  example.  Figure  10  shows  the 
orbits  twith  data)  passing  near  Ft.  Churchill  in  three  consecutive  24-hour 
periods.  A  data  set  was  created  which  contained  all  the  data  available  within 
"boxes",  of  size  six  degrees  latitude  by  bO  degrees  longitude,  centered  on  each 
ot  lb  rocket  stations.  (These  lb  included  the  13  of  Table  1  plus  3  Sll  stations, 
which  had  little  rocket  data  and  were  used  only  to  check  the  time  consistency 
of  SCR  data).  From  this  SCR  time  series  it  was  possible  to  interpolate,  for 
each  rocket  station  and  for  each  day,  a  "daytime"  radiance  from  the  ascending 
orbits  and  "nighttime"  radiance  from  the  descending  orbits,  as  follows.  When 
two  orbits  surround  the  station,  a  station  value  was  interpolated  from  each 
orbit's  average  within  the  box.  When  only  one  orbit  was  available,  the  orbit 
box  average  was  used.  For  each  radiance  channel,  separate  time  series  were 
obtained  from  daytime  and  nighttime  conditions  at  each  rocket  location;  four 
separate  time  series  were  therefore  obtained  for  SCR-A,  and  eight  for  SCR-U. 

The  time  series  of  the  daytime  and  nighttime  radiances  were  plotted, 
along  with  meteorological  parameters  from  the  actual  rocket  observations. 

These  time  plots  were  very  useful  in  identifying  the  many  short  periods  of 
bad  data,  especially  for  SCR-A.  New  procedures  were  then  devised  to  re-clean 
the  data  (these  have  been  discussed  in  Section  IV).  Several  rocket  data  wore 
also  identified  as  erroneous,  and  were  eliminated  from  further  consideration. 

In  our  previous  work  with  SCR-A  data,  we  found  that  large  differences 
sometimes  existed  between  the  daytime  and  nighttime  radiances  at  a  given 
place,  particularly  in  Channel  B.  These  differences  were  larger  than  one 
would  expect  from  diurnal  atmospheric  changes,  and  were  known  to  change  with 


_ _ 


season  anil  latitude.  To  remove  this  effect,  monthly  zonal  means  of  the  diurnal 
differences  for  all  channels  of  both  SCR-A  and  SCR-B  were  computed.  Then, 
interpolating  in  latitude  and  in  time,  we  applied  these  differences  to  each 
ascending  datum,  effectively  making  each  radiance  a  nighttime  radiance. 

B.  Regression 

The  time  series  discussed  in  the  last  section  were  interpolated,  in 
time,  to  the  times  of  rocket  firings.  The  interpolation  was  linear,  and  the 
maximum  time  separation  of  SCR  data  (for  the  purpose  of  interpolation  to  a 
rocket  time  between  them)  was  42  hours,  with  one  exception:  If  the  gap  was 
larger  than  42  hours,  SCR  data  from  an  endpoint  of  the  gap  was  assigned  to  the 
rocket  time  if  the  rocket-SCR  time  difference  was  less  than  6  hours. 

In  order  to  develop  consistent  regression  models,  it  was  necessary 
to  combine  the  rocket  data  into  3  sets  which  represented  latitude  regions  near 
bON,  30N  and  ION,  termed  "high", "mid" ,  and  "low"  (Table  1).  Further,  the  sta¬ 
tion  groups  were  partitioned  into  warm  and  cold  "seasons"  of  generally  six 
months,  which  were  chosen  to  take  into  account  real  atmospheric  changes  while 
allowing  for  changing  instrument  characteristics.  One  of  the  seasonal  parti¬ 
tions  was  chosen  as  17  April  1972,  corresponding  with  the  SCR-A  calibration 
change  mentioned  earlier. 

Several  functions,  X,  of  the  SCR  radiances  were  correlated  with  each 
desired  atmospheric  parameter.  These  SCR  functions  included  the  radiances, 
their  squares,  square  roots,  and  fourth  roots,  as  well  as  the  products  and 
the  ratios  of  the  radiances. 

In  the  regression  development  a  standard  screening  procedure  was 
used.  First,  for  each  atmospheric  parameter  (predictand),  the  most  highly  cor¬ 
related  function,  Xj ,  was  identified  and  the  percent  of  predictand  variance 
which  it  explained  was  calculated.  Then  the  effect  of  Xj  was  removed  from  all 
the  other  correlations  and  the  next  most  highly  correlated  function,  X,,  was 
identified  and  its  contribution  to  additional  reduction  of  predictand  variance 
was  calculated.  The  procedure  was  continued  for  a  third  and  fourth  function 
but  it  was  found,  for  both  SCR-A  and  SCR-B,  that  additional  functions  never 
led  to  a  significant  improvement  over  a  two- function  model.  The  model  was 
therefore  restricted  to  two  functions,  X^  and  X^. 


The  form  of  the  model  used  was 


P  -  P  +  At  .  (Xt  -  Xj)  +  A2  .  (X  -  X  )  (1) 

where  P  is  the  predictand 

P  is  the  mean  value  of  the  predictand 

A.  is  the  coefficient  of  the  i ' th  predictor,  X.  -  X.,  and 
1  i  i 

X^  is  the  mean  Value  of  the  i ’ th  function. 

The  coefficients  A  were  calculated  from: 

r(P,X1)Z  -  r(P,X2)  .  r(X1,X2) 

1  -  r  (X1,X2)2 

r(P,X2)2  -  r(P,X1)  .  r(X1>X2) 

2  o(V  l  -  r(X1,X2)2 

where  o  denotes  standard  deviation  and  r  denotes  correlation  coefficient. 

Such  models  were  generated  for  each  predictand,  for  each  latitude 
region  and  season.  One  problem  was  that  two  functions  which  explained  the  most 
variance  of  a  particular  predictand  were  not  generally  the  same  between  lati¬ 
tude  regions  or  seasons.  However,  in  almost  all  cases,  nearly  as  much  variance 
could  be  explained  by  functions  which  were  judiciously  specified  so  that  they 
varied  smoothly  in  latitude  and  season.  By  varying  smoothly  is  meant  that  one 
of  the  functions  X  must  be  the  same  for  adjacent  latitude  regions  in  the  same 
season  or  at  adjacent  seasons  in  the  same  latitude  region.  (Due  to  the  change 
in  character  of  the  SCR-A  radiances  after  17  April  1972,  smoothness  in  this 
time  transition  was  not  required).  Although  by  specifying  the  model  some  theo¬ 
retically  explainable  variance  was  lost,  this  method  was  preferred  because  it 
insured  smoother  derived  meteorologicai  parameters  across  time  and  latitude 
boundaries.  Very  little  degradation  of  results  for  individual  seasons  or  lati¬ 
tude  regions  was  caused  by  the  adoption  of  this  restriction.  In  the  two  cases 
where  the  penalty  for  imposing  smoothness  was  judged  too  high,  the  restriction 
was  not  enforced. 

It  should  be  noted  that  the  secondary  screening  procedure  was  redone 
when  the  spec i f ied  X^  differed  from  the  truly  most  highly  correlated  function. 


A 

1  "(X^ 
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In  this  way  the  optimum  choice  for  X^  was  assured. 


Appendix  A  contains  details  of  the  models.  The  first  group  of  entries 
for  each  level/latitude/season  (i.e.  "Run  1")  shows  the  functions  used,  predic¬ 
tor  coefficients,  and  the  percent  variance  explained  by  the  model.  (See  p.  A-2 
for  complete  details  of  entries.  Also,  note  that  the  results  for  layer  mean 
temperatures  may  be  viewed  in  terms  of  thickness  by  noting  thaL  1°K  21,  27, 

21,  and  27  m  in  the  10-5,  5-2,  2-1,  and  1-.4  mb  layers,  respectively.) 

Also  listed  in  Appendix  A  are  results  of  tests  of  the  model  which 
were  performed  on  independent  data  in  the  following  manner:  For  each  predic- 
tand/ lati tude-region/season,  a  model  (using  the  specified  functions)  was  devel¬ 
oped  from  only  85  percent  of  the  available  rocket-SCR  pairs.  This  model  was 

« 

then  used  to  compute  the  predictands  for  the  remaining  15  percent  of  the  radi¬ 
ances,  which  were  then  compared  with  the  coincident  rocket  parameters.  The  mean 
error  and  the  standard  deviation  of  the  error  are  given  for  five  such  tests, 
each  of  which  used  a  different,  randomly-chosen,  85  percent/15  percent  combina¬ 
tion  of  dependent/independent  data.  It  is  seen  that  neither  the  variance  ex¬ 
plained  nor  the  coefficients  vary  greatly  as  a  function  of  dependent  data  set, 
and  that  the  standard  deviation  of  the  error  of  the  predictions  is  usually  far 
less  than  the  standard  deviation  of  the  dependent  (rocket)  data. 

These  independent  data  tests  were  also  useful  in  determining  whether 
the  prescribed  and  were  actually  as  good  as  their  "percent  variance  ex¬ 
plained"  advertised  them  to  be.  Often,  a  model  which  was  slightly  inferior  in 
terms  of  percentage  variance  explained,  but  which  was  more  consistent  with  re¬ 
spect  to  its  neighbors  (in  latitude  and  season),  could  be  shown  to  be  of  equiv¬ 
alent  quality  when  applied  to  independent  data.  In  these  cases  the  more  con¬ 
sistent  model  was  the  one  finally  used. 

Averaging  of  the  standard  deviations  of  the  errors  of  the  five  inde¬ 
pendent  data  tests  yields  the  following  accuracies:  At  10  and  5  mb,  the  error 
i9  about  3°K,  generally  increasing  upward  to  about  5.5°K  at  .4  mb.  The  error 
of  the  layer  mean  temperatures  is  generally  about  one-half  degree  less  than 
the  errors  at  the  bordering  levels. 
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As  NMC  LO  mb  temperature  is  available  north  of  20°N,  it  was  realized 
it  should  be  a  good  predictor  of  10-5  mb  mean  temperature  and  other  parameters 
as  well.  However,  since  it  was  not  always  available,  and  never  for  the  tropics 
and  SH,  use  of  NMC  data  was  first  excluded  from  the  regression.  Later  it  was 
felt  better  to  utilize  whatever  data  did  exist  lo  optimize  the  regression,  even 
if  for  only  north  of  20°N.  A  new  set  of  models  was  developed,  then,  using  NMC 
10  mb  temperature  as  a  possible  predictor. 

However,  it  seems  that  NMC  and  MRN  10  mb  temperatures  are  not  always 
well  correlated,  and  after  investigation  it  was  decided  to  use  MRN  as  the  stand¬ 
ard.  Appendix  B  shows  details  of  the  performance  of  the  regressions  when  and 
where  NMC  10  mb  temperature  is  a  useful  predictor;  at  other  levels  and  seasons, 
Appendix  A  results  are  preferable.  This  shows  that  the  prediction  of  tempera¬ 
tures  at  10  mb  and  in  the  10-5  mb  layer  are  improved  by  about  one-fourth  degree 
due  to  use  of  NMC  10  mb  temperatures.  It  is  believed  that  the  accuracies 
presented  in  Appendices  A  and  B  are  quite  good,  especially  considering  the 
errors  present  in  the  data  before  cleaning. 

VI.  APPLICATION  TO  SCR-B  ORBITAL  DATA 

The  models  discussed  in  the  last  section  were  used  to  calculate  atmos¬ 
pheric  parameters  from  SCR-B  orbital  data.  A  full  set  of  atmospheric  parameters 
was  computed  for  every  major  frame  having  data,  with  the  data  first  being  cor¬ 
rected  for  zonal  monthly  mean  diurnal  radiance  differences,  as  explained  in 
Section  V- A.  Several  decisions  were  required  concerning  smoothing,  use  of 
10  mb  (NMC)  data,  and  Southern  Hemisphere  processing,  and  these  are  outlined 
below. 

A.  Smoothing 

It  was  desirable  to  insure  smooth  transitions  of  computed  atmos¬ 
pheric  parameters  across  time  and  latitude  boundaries.  The  time  boundaries 
are  shown  in  Appendix  A;  the  latitude  boundaries  in  the  NH  were  chosen  to  be 
22.5°N  and  A2.5°N  (SH  processing  is  discussed  later).  Smoothness  in  latitude 
was  accomplished  by  1)  computing  the  parameter  using  both  sets  of  statistics 
within  a  +2*  degree  "window"  of  the  boundary,  and  2)  forming  an  appropriately 
weighted  averaged,  based  on  location  relative  to  the  boundary.  Smoothing  across 
time  boundaries  followed  the  same  procedure,  with  the  "window"  being  +A  days. 


B.  Use  of  10  mb  (, NMC )  Data 

Temperature  and  height  at  10  mb  on  the  1977-point  NMC  grid  were  avail¬ 
able,  and,  when  required,  10  mb  temperature  was  used  as  a  predictor  ol  parameters 
north  of  20N.  To  10  mb  height  can  be  added  the  layer  thicknesses  to  determine 
heights  of  5,  2,  1,  and  0.4  mb.  Since  the  NMC  grid  was  only  available  once  a 
day  (12Z),  10  mb  fields  had  to  be  (linearly)  interpolated  in  time,  as  well  as 
space,  to  the  sub- sa tel  1 i te  points.  The  maximum  gap  over  which  interpolation 
was  allowed  was  48  hours  (1  missing  day);  for  longer  gaps  the  regression  models 
not  using  10  mb  fields  were  used. 

C .  Southern  Hemisphere  Processing 

No  extratropical  SH  rocket  stations  had  enough  observations  to  allow 
stable  regressions  with  SCR  data.  It  was  therefore  necessary  to  apply,  in 
the  SH  mid- latitudes  (22.5LS  to  42.5°S)  and  high  latitudes  (poleward  of  42.5°S), 
models  based  on  NH  regressions  in  the  appropriate  time  of  year  (i.e.  winter  or 
summer  half-years).  Of  the  four  NH  regression  periods,  only  the  middle  two 
(16  April  -  15  October  1973  and  16  October  1973  -  15  April  1974)  were  used  for 
this  purpose. 

In  the  SH  tropics,  the  same  regression  model  was  used  as  was  used 
in  the  NH  tropics.  No  6-month  time  shift  was  applied,  since  the  amplitude  of 
the  annual  wave  in  temperature  is  small  (less  than  2°K  from  10  mb  to  1  mb  and 
less  than  4°K  at  .4  mb)  in  the  tropics  (Reference  7). 

Procedures  to  assure  smoothness  over  the  time  and  latitude  boundaries 
were  the  same  for  the  SH  as  for  the  NH. 

VII.  RESULTS 

The  SCR-B  radiances  have  been  used  to  estimate  the  temperature  at  10, 

5,  2,  1,  and  0.4  mb  and  the  mean  temperature  and  thickness  of  the  layers  be¬ 
tween  these  surfaces,  along  the  orbit,  for  the  Northern  and  Southern  Hemispheres. 

In  the  next  phase  of  this  work,  the  heights  of  the  same  surfaces  at  lati¬ 
tude-longitude  grid  points  will  be  computed  for  use  in  obtaining  planetary  waves, 
for  the  five-year  period  of  record  of  SCR-A  and  SCR-B.  This  will  require  care¬ 
ful  space-and-time-interpolation  of  all  the  data. 
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Table  l 


Rocket sonde  Stations  Used  tor  Regressions  with  SCR-A  and  SCk-B 


Latitude 

Grouping 

Station 

No. 

Station  Name 

Lat i tude 

Longitude 

04202 

Thule 

76. 6N 

68. 8W 

/ 70266 

Ft.  Greely 

64.0 

145.7  \ 

High 

\ 70192 

l'oker  Flats 

65.0 

147.5  / 

72913 

Ft.  Churchill 

58.7 

93.8 

72124 

Primrose  Lake 

54.8 

110.1 

72402 

# 

Wallops  Is. 

37.8 

75.5 

72391 

Pt.  Mugu 

34.1 

119.1 

Mid 

72269 

White  Sands 

32.4 

106.5 

74794 

Cape  Canaveral 

28.5 

80.5 

91162 

Barking  Sands 

22.0 

159.8 

78861 

Antigua 

17.2 

61.8 

/  78783 

Ft.  Sherman 

9.3 

80.0  \ 

Low 

\ 78801 

Ft.  Sherman 

) 

91366 

Kwa  j  a  1 e i n 

8.7 

167. 7E 

61902 

Ascension 

8.  OS 

14. 4W 
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WtICHTING  FUNCTIONS  FOR  NIMBUS  4 


vriCfiriNU  FUNCTIONS  TOR  NIMBUS  S 


Figure  1.  Weighting  functions  for  (a)  the  Nimbus  4  SCR  and  (b)  the 
Nimbus  5  SCR.  (From  reference  6). 
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r lgure  6.  (a)  Imbalance  corrected  channel  A  radiances  from  part  of 
Nimbus  4  readout  orbit  1479. 
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Figure  6.  (b)  Radiances  of  6a  smoothed 

twice  with  a  five-point  running  mean. 


Figure  6.  (c)  Radluncea  of  6a 

smoothed  four  times  with  a  five-point 
running  mean. 
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Figure  8.  Examples  of  filtering. 
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APPENDICES 


Statistical  accuracy  of  temperature  retrievals  from  SCR-B  radiances 
at  meteorological  rocket  network  locations,  not  using  10  mb  NMC 
temperatures.  These  show  the  regressions  used  for  the  regions 
south  of  20°N. 

Same  as  A,  for  those  atmospheric  parameters,  latitudes,  levels,  and 
periods  for  which  10  mb  temperatures  improved  the  retrievals,  i.e. 
where  "2"  is  indicated  as  a  predictor.  These  regressions  were  used 
north  of  20°N. 


EXP1.ANAT0RY  notes  for  appendices  a  and  b 


Accuracy  of  Temperat  u it  Ret  r  i  c vals  from  St~K-»  Radiances  at  M RN  Locations 


1.  Regross  ton  models  for  temperature  and  mean  temperatures,  developed  using  all 
available  rocket  data,  are  shown  in  run  1.  These  models  were  then  tested  as 
fol lows : 

a.  For  each  station  group  and  period,  a  randomly  chosen  15  percent  of  the 
available  MRN  observations  were  set  aside  to  serve  as  an  independent  test 
set,  and  the  remaining  85  percent  were  used  to  develop  a  regression  model. 

b.  This  regression  model  was  applied  to  the  radiances  occuring  at  the  other 
15  percent  of  the  MRN  locations  and  times,  with  a  resulting  mean  error 
and  standard  deviation  of  the  error  as  shown  in  the  listing. 

C.  Steps  a  and  b  were  repeated  for  four  more,  different ,  sets  of  random 
data.  These  five  verification  tests  are  shown  under  runs  2-6. 

2 .  Explanation  of  printout  tables: 

A  B  D  F  I 

C  E  G  J 

H  K 


A  «  Run  number 

B  »  Primary  predictor  (see  3  below)  used  in  the  model 
C  *  Secondary  predictor  used  in  the  model 

D  *  Variance  explained  by  a  model  which  uses  only  the  primary  predictor 

E  »  Variance  explained  by  the  full,  2  -  predictor  model 

F  «  Mean  of  the  rocket  data  used  in  the  model,  P  (see  eq.  (1)  in  text) 

G  *  Coefficient  of  the  primary  predictor, 

H  =  Coefficient  of  the  secondary  predictor,  A9 

I  *  Number  of  observations  (run  l),  number  of  independent  cases  tested 
(runs  2-6) 

J  =  Mean  error  of  the  independent  cases  (  K).  Applies  only  to  runs  2-6 
K  =  Standard  deviation  (°K)  of  all  the  rocket  data  (for  run  1);  standard 
deviation  (°K)  of  the  error  for  the  independent  test  cases  (for  runs 
2-6) 


3.  Explanation  of  predictors.  (Note:  Rj  =  ch  B12  radiance,  R-,  =  ch  B23, 
R3  -  ch  B34 ,  R4  =  ch  B4.) 


NMC  10  mb  temp 


3  ■  R-  .  R. 
3  4 


4  -  Rl/R3 

5  -  Rj/R4 


4.  Station 


6 
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10  -  R, 


11 


12  =  R, 


13 


R, 


groups  on  each  page  are  for  high  (60  -75  N) ,  middle  (30  -40  N),  and 
low  (10°S-20  N)  latitudes,  arranged  from  top  to  bottom.  From  left  to  right, 
are  6-month  periods  (year/day).  Temperatures  for  5  levels  are  shown  first, 
followed  by  4  mean- layer- temperatures . 
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